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Abstract. An outline of theoretical estimates is given regarding the dependence of the 
value of the diffusion coefficient D on the size R of the remainder of the normal form in 
doubly or simply resonant domains of the action space of 3dof Hamiltonian systems. 

1. INTRODUCTION 

In the works of Froeschle et al. (2000, 2005, 2006), Lega et al. (2003), and Guzzo et 
al. (2005) precise numerical estimates were given of the critical threshold (in terms of 
a small parameter e) at which one has the onset of the so-called 'Nekhoroshev regime' 
(Nekhoroshev 1977, Benettin et al. 1985, Lochak 1992, Poshel 1993) in conservative 
systems of three degrees of freedom. In the same works the local speed of Arnold 
diffusion, as well as the exponents appearing in the associated laws of the Nekhoroshev 
theory, were estimated. Other studies by the same group are Guzzo et al. (2006, 
effects of non-convexity) and Todorovic et al.(2008, diffusion in 'a priori unstable' 
systems). Earlier studies are: Kaneko and Konishi (1989), Wood et al. (1990), Dumas 
and Laskar (1993), Laskar (1993), Skokos et al. (1997), Giordano and Cincotta (2004). 

In a recent study (Efthymiopoulos 2008), we used a computer program to carry on 
Hamiltonian normalization up to a sufficiently high order, and found that in a case of 
simple resonance the size R of the optimal remainder of the normalized Hamiltonian 
(which turned to be exponentially small in l/e a , a ~ 0.25) scales with the diffusion 
coefficient of numerical experiments as D a R 3 . This scaling law was presented as 
an empirical fact, no theory behind being suggested. Furthermore, no investigation 
was made of what happens in cases of double resonance. The above two questions are 
briefly addressed in the sequel. In particular, theoretical estimates of the scaling law 
D(R) are outlined, first in the case of double and then of simple resonances. Details 
are deferred to a paper under preparation. 

2. OUTLINE OF ESTIMATES ON THE LAW D(R) 
We shall refer to Hamiltonian systems of three degrees of freedom 

H(I,<t>) = H (I) + eH 1 (I,<t>) (1) 
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where I = (J±, I 2l I^), 4> = {<fii,4 > 2, < p2), Hq satisfies appropriate non-degeneracy and 
convexity conditions, and H is analytic in a complexified domain of the actions and the 
angles. The frequencies are oj{I) = WiHq(I). Owing to the analyticity condition, the 
Fourier series H\ = hk{I) exp(ifc • 0) yields coefficients hk decaying exponentially 
with the modulus of the vector fc, namely the bound |/ifc(/)| < Aexp(— \k\a) holds with 
\k\ = | + \k 2 \ + l^l and A, a positive constants. Consider an open domain W/,,d 
in the action space, of size D, centered around some central value /*. Expanding Hq 
around i* leads to 

H = #o* + oJ* ■ J + -M^JiJj + ... (2) 

where J = I — /*, a;* = V/i?o(**), My* are the entries of the Hessian matrix of Hq 
at /* (which satisfy restrictions imposed by the convexity condition) . Assuming some 
truncation order K c in Fourier space, we consider the case in which there are two 
linearly independent integer vectors k^\ k^ both satisfying the conditions |fe^| < 
K c , feW ■ u>* = 0, j = 1,2. A normalization of the Hamiltonian JT]) valid in Wd i„ 
results in that all the resonant trigonometric terms, of the form exp(z((<?ifeW +q 2 k^)- 
4>), with qi,q 2 integers, \q\k^ +q2k^)\ < K c , survive in the normal form. Choosing 
a vector m satisfying m ■ k^ = 0, j = 1, 2, and setting J; = k^I^ + fc^ jjj + rriilp, 
4>ii = ' 4>i 4>f = n% ■ (f>, i = 1, 2, 3, j = 1, 2, the normalized Hamiltonian takes the 
form (in new canonical variables) 

H norm = Z R (l R X) , 4 2) , I F , 0« , <f$ ) + R (/« , if , I F , 0g) , 4$ , 4> F ) . (3) 

Two facts are relevant about ([3]): i) the angle <pp is ignorable in the normal form Zr. 
Therefore If is an integral of the Hamiltonian flow under Zr alone. It follows that, for 
different label values of If, Zr alone defines the dynamics of a system of two degrees 
of freedom, ii) The (optimal) remainder is exponentially small R = 0(exp(— 1/e)), 
thus it only slightly alters the dynamics due to the normal form. 

Now, the convexity condition ensures that there is a linear canonical transforma- 
tion (4\J%\ 

'> 4>%' ') — * J'li 02) sucn tnat m tne new variables the Hamil- 
tonian reads Z R = \ [Ai(J[ - c FA ) 2 + A 2 (J 2 - cp 2 ) 2 ] + ... + eJ2 n C„{J[,J 2 ;e,I F ) 
x exp(in • 0'), the constants Ai, A 2 being positive and depending on My*, and 
c_f,i,cf,2 depending on My* as well as on the value of the integral I F . Apart from 
trivial modifications, the dynamics of Zr is then the same as under the simplified 
model 

Z R = i [(J[ - c F1 ) 2 + (J* - c F2 ) 2 ] +eJ2C n (Ji J 2 ; e, I F ) exp(m • 0') (4) 

n 

which we now consider in some detail. 

For any fixed value of the angles (^ ,4>' 2 ), the constant energy condition Zr = E 
implies that the actions J[,J 2 he on a circle of radius \J2E centered at (J[,J 2 ) = 
(cf,1j cf 2 )- This is shown schematically in Figure 1. As the angles <j>' change value, the 
trigonometric terms in Q can only induce a O(e) change in the radius of the circle, 
which we shall temporarily ignore. On the other hand, the lines vn! x J' x + m' 2 J 2 = 
0, mf ■ n = 0, where n is the integer vector of any trigonometric term appearing 
in define the ID resonant manifolds of the Hamiltonian ([!]) in the reduced 2D 
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action space. By construction, the Hamiltonian @ contains only a finite number of 
harmonics with coefficients C n scaling as |C n | ~ e exp(— cr|n|). It follows that there 
is a finite set of resonant lines in Fig.l, all passing from the center of the circles of 
constant energy. Three such lines are shown schematically (black solid) . The pairs of 
dashed lines defining the boundaries of the zones around the resonant lines correspond 
to the limits of the separatrix width of each resonance, which is of the order 

AJ n ~(eexp(-a|n|)) 1/2 . (5) 

The smaller the value of |n| the larger the width of the resonance. Thus, in Fig.l we 
have |m| < \n 2 \ < \n 3 \. 

The resonances are well separated if the energy E is large enough (big circle in 
Fig.l, at E — E\). Since the entire set of resonances cover an 0(e 1 / 2 ) length of an 
arc on a circle of constant energies, significant resonance overlap occurs if E = 0(e) 
or smaller (small circle for E — E 2 in Fig.l). 

So far we have neglected the effect of the remainder. This is manifested by causing 
a slow drift of the normal form energy value Enf of the orbits run under the full 
Hamiltonian, compensating the drift in the remainder value so that the total energy 
E = Enr + R be a constant. In Fig.l this drift is shown by curly curves. An orbit 
found at a given time near the separatrix of, say, resonance 1, with Enf = E±, has 
a certain probability to slowly drift in the action space outwards (Enf increases) or 
inwards (Enf decreases). Assume the latter case. Then, after a long time the orbit 
will be found touching the inner circle E = E 2 , where all resonances overlap. If, still 
later, the drift is reversed (again with a certain probability), Enf will increase and 
the orbit will return to the circle Enf = E\, but this time lying, in general, on the 
separatrix layer of a different resonance, say resonance 2. 

This qualitative picture of the diffusion in resonance junctions should, of course, be 
substantiated by calculating the invariant manifolds of the 2D tori lying at the borders 
of the resonances of the full Hamiltonian system. If there are heteroclinic intersections 
between these manifolds, an orbit started on one manifold is guaranteed to undergo 
the type of drift motion described above. However, verifying this fact numerically 
poses a formidable task, equivalent to probing the very mechanism of Arnold diffusion. 
This restriction notwithstanding, we proceed in a plausible quantitative estimate of 
the speed of the diffusion on the basis of numerical experiments having observed, 
instead, the drift in the action space directly (as e.g. in Lega et al. 2003). The 
indication from such experiments is that the drift in the action space can be modeled, 
at least locally, like a normal diffusion process. Assuming this true, one may ask how 
long it will take for, say, the curly itinerary of Fig.l to be realized. Attributing the 
drift to the remainder, the per step change of the Enf value of an orbit crossing 
the junction is dE ~ R. By Fick's law, one then has after N steps a total change 
AEnf ~ N x / 2 dE — N X / 2 R. By © the radius of the outermost circle touching the 
junction's limits is A J = 0(e 1 ^ 2 ), implying that the total drift in the value of Enf 
is AEnf = 0(e). Putting these relations together one then finds for the diffusion 
coefficient D ~ AJ 2 /N the relation D ~ e^R 2 . Note that, since R = 0(exp(-l/e)) 
one has, by de l'Hopital's rule that D ~ R 2 — + in the limit e — > 0. 

A similar calculation can be made in the case of simple resonances, if we take into 
account that the simply-resonant domains are, in fact, composed also by sub-domains 
VVd.i, crossed by multiple resonances, the difference from doubly resonant domains 
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Figure 1: The circles of constant normal form energy E^f and the drift of the orbits 
in resonance junctions (see text) 

being essentially that only one of these resonances satisfies \k^\ < K c , while all 
the others satisfy lA;*- 2 ^ > K c . Since all apart one resonances are contained in the 
remainder function, it follows that the total travel to cross a resonant junction in 
the action space is now of order i?, 1 / 2 (instead of e 1 / 2 in double resonances). Setting 
A J 2 ~ R we then readily find D ~ e~ 2 R 3 . 

In conclusion, we predict that in systems like Jl]) the local value of the diffusion 
coefficient D scales with the size of the optimal remainder function R of the local 
normal form as a power law D ~ R b , the estimate b « 2 holding close to doubly- 
resonant domains and b m 3 close to simply-resonant domains. It would be of interest 
to check these theoretical estimates against detailed numerical experiments. 
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